Anomalous Negative Magnetoresistance Caused by Non-Markovian Effects 
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A theory of recently discovered anomalous low-field magnetoresistance is developed for the system 
of two-dimensional electrons scattered by hard disks of radius a, randomly distributed with concen- 
tration n. For small magnetic fields the magentoresistance is found to be parabolic and inversely 
proportional to the gas parameter, 8p xx /p ~ ~{uj c t) 2 /na? . With increasing field the magnetore- 



sistance becomes linear 5p xx /p 
simulations. 



-lu c t in a good agreement with the experiment and numerical 
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It is well known that in the Boltzmann-Drude ap- 
proach the longitudinal resistivity p xx of a degenerate 
two-dimensional (2D) electron gas does not depend on 
the transverse magnetic field B. Therefore, the known 
mechanisms of magnetoresistance (MR) involve either 
quantum interference effects or classical non-Markovian 
memory effects, which are not captured in the Boltz- 
mann picture. The MR, arising from quantum effects 
was discussed in a great number of works (see for re- 
view Ref. 0). The role of classical memory effects was 
underappreciated for a long time, though several theoret- 
ical works pointed out at the importance of such effects 
for magnetotransport 0,1110]. The interest to the prob- 
lem of classical MR has sharply increased in recent years, 
starting with Ref. [f| , where it was shown that effects of 
"classical localization" may lead to the exponential sup- 
pression of electron diffusion at large B. This work was 
followed by a series of works 0, Hll El El El El , 
discussing different aspects of classical 2D magnetotrans- 
port in strong magnetic fields. 

In this paper we focus on a mechanism of low field clas- 
sical MR specific for systems of strong scatterers. This 
mechanism is connected to the memory effects due to 
backscattering events. The corresponding corrections to 
the conductivity are small in the parameter a/ 1, where a 
is the characteristic size of the scatterers and I is the mean 
free path. Nevertheless, the dependance of these correc- 
tions on the magnetic field turns out to be very sharp, 
resulting in the MR anomaly. The anomaly was discov- 
ered in recent numerical simulations 15] where the MR 
in a system of 2D electrons scattering on randomly dis- 
tributed hard disks was studied. This system is usually 
referred to as the Lorenz gas and is the simplest model 
of the 2D electron gas with strong scatterers. In the fol- 
lowing we restrict our considerations to this model. The 
generalization of our results to other models of strong dis- 
order is straightforward. Magnetotransport in the Lorenz 
gas is characterized by two dimensionless parameters: 
= lo c t, and the gas parameter 0o — a/ 1 = 2na 2 . Here 
a is the disk radius, n is disks' concentration, oj c is the 
cyclotron frequency, r = l/vp is the mean free time and 



I = X/2na is the mean free path. The anomaly was ob- 
served in the case <C 1, /3q <C 1. Both the numerical 
simulations and the qualitative analysis of 1^ indicated 
that at zero temperature, T, the MR can be expressed in 
terms of a dimensionless function f(z) via 
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where p is the resistivity for B = 0. Numerical re- 
sults 15'] suggest that /(z) ~ z as z ^ 0, yielding 
5 pxx I P ~ — |w c |t. The latter expression is in a very good 
agreement with experimental measurements of MR in a 
random antidot arrays |l6j . It is anomalous in two senses. 
First, it has a non-analytic dependence on the magnetic 
field. Second, it does not vanish in the limit of vanishing 
00, which is normally regarded as the expansion param- 
eter for the corrections to the Drude-Boltzmann picture. 
This intriguing behavior calls for a rigorous analytical 
theory of the effect, which would establish Eq. Q]and en- 
able one to derive the analytical expression for function 
/. In this letter we present a theory of the anomaly and 
give an expression for f(z). We find that for z < 1, f{z) 
is linear in agreement with numerical experiment, but at 
very small z < 0.05 crosses over to a quadratic depen- 
dence. Thus, for — > 0, Eq. ^yields 5p xx /p ~ ~0 2 /0 o . 
The limit 0o — * should be taken with care. While the 
small expansion seems to be singular as a function of 
0o, the region of where this expansion is valid shrinks 
as 0o — > 0. For z — > oo, / saturates at some constant 
value. Therefore, the full variation of 5p xx /p is of the 
order 0q. In other words, the anomalous MR is strong 
but it exists in a small region of magnetic fields. 

In E3 a mechanism of MR connected with mem- 
ory effects arising in backscattering events was pro- 
posed. It has a close relation to the well known non- 
analyticity of the virial expansion of transport coefficients 
E3 El El El El, which we briefly recall. For B = 
the leading nonanalytic correction to resistivity, 5p, is 
due to the processes of return to a scatterer after a sin- 
gle collision with another scatterer [Fig. 1(a)]. The rela- 
tive correction, Sp/ p, is proportional to the correspond- 
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FIG. 1: Backscattering process responsible for leading non- 
analytic contribution to the resistivity at B = (a). For 
B / 0, the overlap area, Sb, between two corridors is small 
at large B (b). For <j> = 0, Sb decreases with B (c). For 
4> 7^ and small B the values of Sb ~ So for time reversed 
trajectories have opposite signs (d,e). 



ing backscattering probability, given by the product of 
e~ r / l d(j)dr/l (which is the probability to reach scatterer 
2 without collision and scatter in the angle d(f>) and the 
probability p to return without collisions from 2 to 1 (here 
I is the mean free path). Assuming p = exp(— r/l) and 
integrating over intervals < <f> < a/r, a < r < oo, one 
obtains EMI 13 IH E2 Sp/p ~ A) MV 2 A))- 

In Ref. [l5j it was shown that the probability p is ac- 
tually larger than exp(— r/l) because the passage of a 
particle from 1 to 2 ensures the existence of a corridor of 
width 2a free of the centers of the disks. This reduces the 
scattering probability on the way back, yielding p(r, (f>) — 
exp(— r/l + nSo(r, </>)), where So(r, (j)) = 2ar — r 2 |</>|/2 is 
the area of the overlap of the two corridors [Fig. 1(a)]. 
For example, for <fi = 0, we have So = 2ar and p = 1. 
Physically, this means that the particle is unable to scat- 
ter, since it travels back along the same path. Taking 
into account the effect of "empty corridor" , we get 
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where C is a constant of the order of unity. Thus, for 
B = this effect simply changes the constant in the 
argument of the logarithm. 

The key idea suggested in 0] was that for B ^ 
the area of the overlap of the two corridors, Sb, sharply 
depends on B, resulting in the observed MR. Indeed, it 
is seen from Fig. 1(b) that for (3 > f3o, Sb — > 0, resulting 
in sharp negative MR 
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The following qualitative explanation of the observed lin- 
ear MR was presented in Ref. 0]. The value ti(Sb — So) 
was estimated for <j> — [Fig. 1(c)] to the first order in 
B as —nr 3 /R c = —r 3 /2alR c , where R c is the cyclotron 



radius. Assuming that this estimate also works at <f> ^ 
and expanding e nSs — e nS ° to the first order in B, one 
gets 6p xx /p ~ -l/R c = -u c t. 

In fact, the physical picture of the phenomenon is more 
subtle. The contribution of any trajectory with <j) ^ is 
cancelled in the first order in B by the contribution of the 
time-reversed trajectory, since the values of Sb — So are 
opposite for these paths [ Figl. (d),(e)]. The cancellation 
does not occur only at very small (j> ~ /3. The integration 
in Eq. Elover (j) < (3 yields 5p xx /p ~ —f] 2 /f] . Larger 
values of <fi also give a quadratic in (3 contribution to 
the MR. This contribution is positive and comes from 
the second order term in the expansion of e nSs — e nS ° 
in B. It follows from our results [Eqs. Q,©] that the 
contribution of small angles is dominant resulting in a 
negative parabolic MR [22j. We find that the parabolic 
MR crosses over to linear at very small (3 m 0.05/9o, which 
explains why the parabolic MR was not seen in numerical 
simulations 0] and experiment [lif . 

Next we sketch our calculations. We consider the 
Lorentz gas at T = 0, assuming that |3<1, (3 <C 1. 
In this case j^, Sp xx /p = —{D — D )/D Q , where D = 
i>i?Ztr/2 is the Drude diffusion coefficient for B = 0, 
l iT = 31/A = 3/8na is the momentum relaxation length 
and D is given by 
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\ «<v(0)v(t)) = \ J drdv(G)vv . (4) 



Here G = G(v, vo, r) is the Fourier transform (at u = 0) 
of the retarded Green's function of the Liouville equation 
and (...) stands for the averaging over the positions of 
the disks. The equation for G reads 



(Zo-T_-T + )G = 5(r)<5(v-v ), 



(5) 



where Lq = vd/dr— w c [vx9/9v] is the Liouville operator 
of the free motion in the magnetic field. The interaction 
with disks is written in Eq. in the form of a collision 
integral |18l \24\ . The scattering operators transform 
arbitrary function /(r,v) as follows, 

f+/(r,v)=p F / dv'a{ip)5(e-e')n + f{r,v% 



f~f(r, v) = -p P / dv'*(ip)S{e - e')n-f(r, v), (6) 



where ri ± = X)i^( r ^ Ri i a)- Here Ri are the po- 
sitions of the disks, delta-function S(e — e') provides 
the energy conservation, pp is the Fermi momentum, 
a(ip) — {a/2)\ sin((/?/2)| is the differential cross-section 
of one disk and (p is the angle between v' and v. The 
vector a = a(v', v) = a(v' — v)/y^2(v 2 — v'v) is pointing 
from the center of a disk to the scattering point at the 
disk surface [Fig. 2). Physically, operator T + describes 
influx of particles with velocity v at the point — a, 
while operator T~ describes the outflux of particles with 



FIG. 2: Scattering of a particle on a hard disk. 
N =1 N = 2,3,... 

' " " \ ' C - ^ \ 

7 \ 7 ' _ \ \ 

vAAA^ X <VWV JVJVL-I- S X N - --— "WW 

± ± +----+ 

(a) ~ " (b) 

FIG. 3: Diagrams, corresponding to the process shown in 
Fig. 1(a). Diagram (a) does not take into account effect of 
"empty corridor" and should be renormalized by (b). 



velocity v at the point Ri + a. The Boltzmann equation 
is obtained from Eq. Elby averaging the Liouville opera- 
tor over the positions of the disks, yielding (T + )/(r, v) = 
p F n J dv'<T(v,v')6(e - e')/(r,v'), (f~) = -1/r. Here 
n = (n^) is the concentration of the disks and 1/r = 
nv2a is the inverse full scattering time. Introducing now 
ST^ = T ± — (T ,± ) and writing a formal solution of J5J, 
G=(6 + L -f- - T+)- 1 as a series in ST*, we get 

(G) ~ Gq + ^ Go(6T a GQ5T J )Go + ■ ■ ■ , (7) 

a, 7 — ± 

where G = (<5 + £o + V T - is the Green's func- 

tion of the Boltzmann operator (here S — > 0). Eq. [7| 
gives a regular way to calculate correlations, which are 
absent in the Boltzmann picture. 

Consider first the case B = 0. Substituting the first 
term in the right hand side of Eq. into Eq. 
we get D = D - The second term in Eq. [7| describes 
the memory-effect due to diffusive returns. As discussed 
above, the main contribution comes from returns after 
a single scattering. This process is described by the di- 
agram Fig. 3(a). The dashed line corresponds to the 
pairings (5T a 5T 1 ) (a, 7 = ±), external wavy lines to the 
diffusion propagators Go- The internal line corresponds 
to the Boltzmann propagator truncated at one scatter- 
ing 6-(f+)G-, where G_ = (L + is the ballis- 
tic propagator and (T + ) stands for one scattering event 
(G_ are shown by solid lines and (T + ) by the cross). 
This diaCTamyields dp/p = -5D/D = (2/3 /3) ln(l/2/3 ) 
H El HI El The terms of the N - th order in 
Eq. contain N pairings (N dashed lines) and are typi- 
cally small as /3q . However, there is a series of diagrams, 
shown in Fig. 3(b), whose contribution is of the order (3q 
19]. The internal dashed lines in this series only contain 
pairings (5T~5T~). Summing the diagrams Fig. 3(b) to- 




FIG. 4: Backscattering process is parameterized by the angles 
ipo,ipf. The magnetic field changes the backscattering angle 
4> = <po + <pf + r/R c - The solid (dashed) line in (a) represents 
electron trajectory for B = (B / 0). Different processes 
contributing to MR are shown in (b)-(e). 

gether with Fig. 3(a), one gets an exact equation 
8p nl t [°°dr _ 2r/l [ 2 * 

P 4t J a r JO JO 

(1 - e lV0 )(l - e W) e nSo(r,0o+0/) = ^ In (jp^} , (8) 

instead of qualitative estimate Eq. [21 Here ipo,iff are 
the scattering angles [Fig. 4(a)], 0o ~ (a/r) cos(y>o/2), 
4>f w (a / r) cos(ip f / 2) and G w 1.8 Thus, addition of 
the series Fig. 3(b) to Fig. 3(a) leads to the following 
renormalization: ln(l/2/3rj) — * ln(G/2/3o)- Physically, the 
series Fig. 3(b) accounts for the effect of the "empty cor- 
ridor" . The N—th order term in this series corresponds 
to N — 1 term in the Taylor expansion of the exp(nSo) in 
Eq. El Four terms in the product (1 - e lVo ){\ - e'^) = 
1 - e lVo - e llp f + e i( - ipa+v i ) correspond to four combina- 
tions of (±, ±) at the ends of external dashed lines in the 
diagrams shown in Fig. 3. They are connected with four 
different types of correlation at a given point r. The dia- 
gram (+, +) [Fig. 4(b)] corresponds to the process, where 
an electron has two real scatterings on a disk placed at 
point r. The diagram (— ,— ) [Fig. 4(c)] does not corre- 
spond to any real scattering at point r. It just allows us 
to calculate correctly the probability for an electron to 
pass twice the region of the size a around point r without 
scattering. To interpret the diagram (+, — ), note that in 
the Boltzmann picture, which neglects correlations, the 
following process is allowed. An electron scatters on a 
disk and later on passes through the region occupied by 
this disk without a scattering [Fig. 4(d)] The diagrams 
(+, — ) correct the Boltzmann result by substracting the 
contribution of such unphysical process. Analogous con- 
sideration is valid for diagram (— , +) shown in Fig. 4[e]. 

For B ^ the sum of diagrams shown in Fig. 3 can 
be expressed as an integral over angles foi^f (scattering 
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angles for B = 0). The only difference from Eq. |HJ is 
that one should replace So — > Sb- For /? <C 1 the overlap 
area can be calculated as Ss(r, <fi) — J Q r dxh(x), where 
h(x) w (2a- |0a;-x 2 /-Rc|) 0(2a- \cj)x - x 2 / R c \) , 9 is 
the Heaviside step function and <f> — 4>q + <f>f + r/R c 
[Fig. 4(a)]. The value of Sp xx /p is obtained from Eq. |H| 
by replacing e nSo to e nS ' B — e nS ° . Introducing dimension- 
less variables T — r/l, z — we get Eq. ^ where 
function f{z) is given by 



3 
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dT 
— f 
T 



-2T 



dtpo / dipf 



, fo + ff \ . > 
cos — - sin 



(f) S in 2 (f)(e--e-). (9) 



Here 



c 



dt 1 



2 

cos((/?o/2) + cos(ifif/2] 



2T 



1 



so 



2 



Function f(z) has the following asymptotics 
/(*) = 



0.33z 2 for z < 0.05 

0.032 (z - 0.04) for 0.05 < z < 2 
0.39-1.3/v^ forz-s-oo. 



(10) 



(11) 



Note that there is a paramctrically small nonanomalous 
correction to Eq. ^due to returns after multiple scat- 
terings, 5p' xx / p i=s —0.2p p 2 0. To compare the results 
of simulations 0] with the theoretical results in a wider 
region of parameters f3, /3q, we substract 5p' xx / 'pfrom the 
numerical curves. Theoretical and numerical | l5| results 
are plotted in Fig. 5. in the universal units, Sp xx /p/3o ver- 
sus z = (3/Pq. It is seen, that the theoretical and numeri- 
cal results are in a ve ry g ood agreement. The comparison 
with the experiment |l6| is more difficult, because of the 
50 % uncertainty in the sizes of the antidots. However, a 
good agreement with the experiment can be achieved by 
appropriate choice of a in the uncertainty interval [l5| . 

Note finally that we fully neglected quantum effects. 
This is possible when a ^> \f\pl (Xp is a Fermi wave- 
length). This criterion ensures that diffraction effects on 
the edges of the disks are not relevant at the scales of the 
order of I. In the opposite case, a <C \/XfI, the diffraction 
should destroy the "corridor effect" , does suppressing the 
anomalous MR. The detailed analysis of quantum effects 
will be presented elsewhere. 

In summary, we have proposed a theory of the negative 
anomalous MR in the Lorenz gas. The analytical expres- 
sion for the MR [Eqs. JTJ), ©, (TTJ] has been derived. 
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FIG. 5: The value of Sp xx /pPo from Eqs. Q, © (solid line) 
shown as a function of /3/(3p together with the results of nu- 
merical simulations presented for different values of /3o 
(triangles for /3o = 0.09 , boxes for /3o = 0.06 , circles for 
/3o = 0.03). Data for all numerical curves are shown for 
/3 < 0.3. Inset: The crossover from quadratic to a linear de- 
pendence at /3//3o ~ 0.05. This crossover was not resolved in 
numerical simulations. 
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